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Abstract — this study explores the capacity of different 
analytic models to predict the vertical velocity profile in 
vegetated flume with submerged vegetation and select the 
suitable model The verification is determined due to a 
comparison between the measured data in vegetated flume 
and the simulated velocity profiles using the Scilab software. 
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I. Introduction 

Vegetation within rivers and flood plains occurs under different 
forms; it can be flexible or rigid and submerged or emergent in 
flows. It plays an important role in the flow structure of many 
streams and rivers and may increase flood risks. Thus, an 
understanding of vegetated flows is necessary to control floods 
and the ecosystem of the stream (Wu and He, 2009; Liu and 
Shen, 2008; Evangelos, 2012). The effects of vegetation on 
flow have been studied for many years by laboratory 
experiments and numerical models with rigid cylinder, flexible 
vegetation prototypes and natural vegetation on open-channel 
flow (Shimizu and Tsujimoto, 1994; Lopez and Garcia, 2001; 
Tsujimoto et al, 1993; Kubrak et al, 2008; Jarvela, 2005; 
Klopstra et al, 1997; Stone and Shen, 2002; Baptist et al., 
2007; Huthoff et al., 2007; Yang et Choi, 2010). Various 
approaches have been proposed are based on two -layer 
approach to describe flow through submerged vegetation. In 
this approach, the flow domain was divided into two layers 
(Figure 1); one through the vegetation called "vegetation layer" 
and the other above it called the "Surface layer". The flow in 
each of the two layers was described separately. The 
logarithmic flow velocity profile is adopted for solving the 
velocity above the vegetation, and the momentum equation 
within the vegetated layer. The continuity of the velocity and 
the shear stress between the two layers is ensured by boundary 
conditions at the interface (Klopstra et al., 1997; Huai et al., 
2009; Jarvela, 2005). For non- vegetated flow, the vertical 
velocity distribution is related directly to the bed shear stress. 
While, for vegetated flow, it is related to the vegetation drag 
because the vegetation roughness is much larger than the river 
bed roughness. Herein, this study compares different analytic 
models based on two-layer approach to analyze the 
performance of these models, for predicting the vertical flow 
velocity profiles through submerged vegetation and to select 
the suitable model to be applied in a real case of river. 
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Figure 1. Velocity profile within and above vegetation 



II. Theoretical Background 

Model 1 klopstra et al.(1997) 

Klopstra et al. (1997), proposed an analytical expression for the 
velocity distribution based on the momentum equation for the 
vegetation layer assuming uniform steady flow. The vertical 
velocity distribution in the vegetated layer is given by the 
following equation: 



u(z) = ^Ciexp(-V2Az) + C 2 exp( A /2Az) + U so 
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Where z is the vertical coordinate, g is the acceleration gravity, 
i is the energy gradient, m is the density of vegetation, D is the 
diameter of plant stems, C D is the drag coefficient, u(z) is the 
flow velocity at the level z, a is a closure parameter derived 
from experimental data and u s0 is the characteristic constant 
flow velocity in non- submerged vegetation. 
The constants Q and C 2 follow from boundary conditions; at 
the bed (z = 0), the bottom shear stress is neglected and the 
flow velocity is assumed to be equal to u s0 , while, at the top of 
the vegetation layer, the shear stress is determined by the 
following expression: 

<h p ) = pg(h -h p )i (4) 

With x is the shear stress, p is the density of water. 

In the surface layer, the velocity follows a logarithmic profile. 

The connection between the boundary conditions at the 
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interface ensures the continuity of the velocity and the shear 
stress between the two layers, and allows the determination of 
the logarithmic law parameters: 



11 * 

u(z) = — ln( 

K 



z-(h D -hs) 



zO 



-) 



(5) 



k is Von Karman's constant (0.41), h s is the distance between 
the vegetation top and the surface layer virtual bed (z 0 <h s <h p ) 
z 0 is the length scale for bed roughness of the surface layer (m) 
and u* is the virtual bed shear . 

Model 2 Righetti and Armanini (2002) 

Righetti and Armanini (2002) developed an analytical two- 
layer model in order to describe a uniform flow condition. Fully 
submerged vegetation is considered as spherical shape. 
Differential equations are obtained by applying the double 
average (in time and space) of the Navier Stockes equations 
neglecting viscosity term. The integration of these equations 
from the free surface flow to the generic level z, leads to the 
following relationships for velocity distribution. 
The vertical variation of the velocity in the vegetation layer is 
expressed as follows: 



h P ,z 



(6) 



u(z) p( ) 
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The only unknown parameter in this expression is the mixing 
length 1. This parameter represents the scale integration process 
of turbulence and the turbulent dispersion. It was determined by 
experimental data. 

In the surface layer, the velocity is defined by the following 
expression: 



u*^ , K(z-(h-hp)) 



u(z) = — ln(l + 

K 



lo 



(7) 



u* is the shear velocity. 

1 0 is the value of the mixing length at the interface level and in 
vegetated layer. 

Model 3 Baptist et al. (2007) 

Baptist et al. (2007) model is based on an analytical solution of 
the momentum balance of flow through and over vegetation. 
The Reynolds stress x is determined using Boussinesq's eddy 
viscosity approach the mixing -length theory to define the eddy 
viscosity: 

Vt (z) =C p lu(z) (8) 

Where v t is the coefficient C p is the turbulent intensity, height - 
averaged over the vegetation height h p and 1 is the mixing 
length. Baptist et al. (2007), proposed an empirical relationship 
for C p and it is expressed as follow: 

C P l = ^(h-h P ) (9) 

The expression of the velocity in the vegetation layer is given 
by the following equation: 



(10) 



u(z) = yu;o + a v exp(y)i 
With: 
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L is the length scale (m) and a v is integration constant. 

For the surface layer, Prandtl's mixing length concept is 

adopted, and the velocity, follows a logarithmic profile: 



u(z) = — ln( ) 



(14) 



K ZO 

Where d is the zero-plane displacement (m), which is located at 

distance from the bed inside the vegetation. 

z 0 is the roughness length and u* is the shear velocity. 

Model 4 Yang and Choi (2010) model 

Yang and Choi (2010), developed a simple relationships based 
on the two -layer approach, and on the momentum equation for 
the vegetation layer assuming uniform steady flow. Yang and 
Choi (2010) neglected the bottom shear stress, assuming that 
the stem drag is the dominant in the vegetation layer. Then, the 
velocity in this case is constant and it is given by the following 
expression: 



Ur 



2ghi 



(15) 



I aCohp 

Where a is the density vegetation (m 1 ). 

In the surface layer, the velocity follows a logarithmic profile 
and depends on the density of vegetation: 



u(z) = ^In(-) + Ui 



(16) 



K Z0 

In this model the roughness length equals to the vegetation 
height (zq = h p ). 

C^ is a constant depend on the vegetation density and it is given 
by the following equations: 

|c„ = l for a^m" 1 ^ 
|c u = 2 for a>5m _1 

u* is the interfacial shear velocity and it's determined by the 
following equation: 

(18) 



u* = A /g(h-h p )i 



III. Results and Discussion 

The verification and the analyzing of these different analytic 
models for predicting the vertical velocity profiles are 
determined due to a comparison between the measured and 
simulated velocities. 

The simulations were carried out by the Scilab software to 
reproduce the vertical profiles of measured velocities in a given 
section. 

Table 1 shows, the experimental data flume available in the 
literature, concerning the free surface flow in presence of rigid 
and flexible, artificial and natural vegetation. 
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Experiments data of Shimizu and Tsujimoto (1994) and Lopez 
and Garcia (2001) were used to verify the capacity of these 
models in the case of rigid vegetation with high and low 
density. Experiments data of Tsujimto et al.(1993) and Jarvela 
(2005) were used to verify the capacity of these models in the 
case of flexible vegetation ( artificial and natural). 
Table 1. Experiments data used for verification of the 
capacity of the three models in prediction of the velocity 
profile 



Author(s) 


Ty 
pe 


h( 
m) 


(m 

) 


D(m 

) 


m 
(m 

2 ) 


C 

D 




Shimizu and 
Tsujimoto (1994) 


R 


0.0 

936 


0. 
04 
6 


0.00 
15 


250 
0 


1 


0.02 


Lopez and Garcia 
(2001) 


R 


0.3 
35 


0. 
12 


0.00 
64 


170 


1. 
1 

3 


0.00 
36 


Tsujimoto et al. ( 
1993) 


F 


0.1 
1 


0. 
06 

2 


0.00 
062 


100 
00 


2 


0.01 


Jarvela (2005) 


N 


0.0 
936 


0. 
04 
6 


0.00 
15 


250 
0 


1 


0.02 



R-Rigid vegetation, F- Flexible vegetation 

The comparison between the measured and calculated velocity 

profiles was illustrated in the following figures: 
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Figure. 2. Measured (Shimizu and Tsujimoto, 1994) and 
simulated flow velocities through rigid vegetation with high 
density. 
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Figure. 3. Measured (Lopez and Garcia, 2001) and simulated 
flow velocities through rigid vegetation with low density. 
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Figure 2 shows simulation results in the case of rigid vegetation 
with high density. In this case, both klopstra et al. (1997) model 
and Baptist et al. (2007) model indicates a performance in the 
prediction of the vertical velocity profile. In contrast, the 
models of Yang and Choi (2010) and Righetti and Armanini 
(2002) demonstrate an under estimation of the velocity. Yang 
and Choi model, assumed that the bed roughness height z 0 
equals to the vegetation height, while, It was supposed smaller 
than the vegetation height in the models of klopstra et al. 
(1997) and Baptist et al. (2007), and this could explain the 
underestimation observed in the surface layer. For Righetti and 
Armanini (2002) model, the difference between the measured 
profiles and the calculated is due to the mixing length 1 0 . This 
parameter was determined by an empirical formula and 
validated only for a limited condition of a laboratory 
experiment. In the case of rigid vegetation with low density, the 
comparison between the simulated and measured profile 
indicates an inefficiency of whole models in the prediction of 
the vertical velocity (Figure 3). These analytic models neglect 
the effect of the bed roughness behind the vegetation roughness 
and this could explain the deviation observed between the 
measured and calculated profiles in this case because Bed 
roughness becomes more important, for sparse vegetation. 
All of these analytical models are set for rigid vegetation and 
it's questioned about their performance in calculating the 
behavior of flexible vegetation. 

Figure 4 and 5 show the comparison between the measured and 
the calculated profiles through flexible vegetation. 
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Figure 4. Measured (Tsujimoto et al., 1993) and simulated 
flow velocities through flexible vegetation. 

In the case of flexible and artificial vegetation (Figure 4), 
calculated profiles by Klopstra et al. (1997) model and Baptist 
et al (2007) model were in good agreement with the measured 
profiles 

The prediction of the velocity profile by these two models is 
also useful in the case of flexible vegetation. However, the 
simulations by Yang and Choi (2010) model, indicate an under 
estimation spatially in the vegetation layer and an over 
estimation was observed with the simulation by Righetti and 
Armanini (2002) model. 

Figure 5 shows the performance of Baptist et al.(2007) model, 
for predicting the vertical velocity in the case of natural 
vegetation, only near the water surface, it is slightly over 
predicting the measured velocity profile. 
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Figure 5. Measured (Jarvela, 2005) and simulated flow 
velocities through natural vegetation 

For flexible vegetation, these models are less accurate than in 
the case of rigid vegetation. The flexibility of the vegetation is 
not even taken into account by these models. This could explain 
the deviations between the predicted and calculated velocity 
profiles. 

IV. Conclusion 

The results of these simulations for submerged vegetation show 
a performance of the two models of klopstra et al. (1997) and 
Baptist et al. (2007) in the prediction of the velocity profiles for 
rigid, flexible and natural vegetation with a high density. In 
contrast, the comparison between the measured and calculated 
profiles by Righetti and Armanini (2002) model and Yang and 
Choi (2010) model, indicates a discrepancy. The use of the 
proposal formula of the mixing length by Righetti and 
Armanini (2002) and the vegetation height as a roughness bed 
height could explain this offset deviation between the measured 
and calculated profiles. In perspective, we will consider tests 
for other models using more experimental data, to verify the 
prediction of the vertical velocity profiles and the mean 
velocity and select most suitable model that will be applied for 
a real cases, the river of Isere in France and Medjerda in 
Tunisia. 
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